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ABSTRACT 

The spcctinim of ocean dynamics lends Itself to convenient subdivision into two 
components in the context of satellite remote sensing. The first is the quasi- 
stationary constituent, while all features which vary with time during the period 
of data acquisition, comprise the second. The precision achieved in each case, 
must be at least d=10 cm through wavelengths of interest, 

Data collected at the Earth's surface — gravity, sea surface heights — cannot 
play a role in the determination of the global gravity field to better than ±0. 3 
mGal (or ±lm) unless assumptions are made about the global characteristics 
of the sea surface topography (SST) — a proposition which is untenable. 

Satellite-determined gravity fields are the only source of data on the gravity 
field which is potentially uncontaminated by data referenced to the ocean sur- 
face. However, the resolution obtained is correlated with the noise level of the 
tracicing and through wavelengths which are a function of spacecraft height, 

Tlie major problems to be overcome at the present time are the improvement 
of the precision of the GEM models to the desired levels and extending the reso- 
lution of the models to shorter wavelengths of practical significance for ocean 
dynamic modelling (at least 5001cm). An appropriately configured GRAVSAT is 
one possible means for obtaining the necessary information with wavelengths 
between 5001cm and 20001cm, 

The definition of time variations in SST is needed for the synoptic modelling of 
ocean dynamics. A Icnowledge of the gravity field is required in this context 
only for recovering the radial component of altimeter spacecraft position globally 
to at least -x-lO cm as the non-tidal variatioii of geoid heights with time does not 
exceed X5 cm, Tlie gravity field model re^juirements in this case are less exact- 
ing in that constituents with wavelengths shorter than those affecting radial com- 
ponents of orbital position with magnitudes less than the desired precision, need 
not be lcno^vn. Tlie effect of the permanent Earth tide on determinations of die 
quasi-stationary components of the spectrum are also examined. 


'•'On leave of absence from the University of New South Wales, Sydney, Australia, 
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THE EARTH’S GRAVITY FIELD AND OCEAN DYNAMICS 


1. GRAVITY ANOMALIES AND THE GRAVITY FIELD 

In the predictable pursuit for Intellectual elegance, physical geodesists have 
sought to formulate solutions of the geodetic boundary value problem In terms 
of surface integrals. Practical considerations have constrained the use of such 
Integrals in numerical evaluations. Prior to 1957, it was hoped to solve these 
Integrals from surface gravity data (e.g. , Molodenskli, et al, , 19G2; HeisUanen 
and Moritz 19G7), No serious attempt was made to assess the level of precision 
attainable oven if a global coverage of data were available. Except for a hand- 
ful of absolute determinations, most gravity values (g) are established using 
differential determina cions. The resulting surface gravity values can be com- 
piled in relation to a world-wide gravity standardization netw'ork {IGSN 71) which 
was established in 1971 with an estimated precision of ±0,2mGal (Notej 1 Gal - 
1cm s'-), as reported in (Morelll, etal,, 1971). The internal precision of 
modern regional gravity networks is likely to be of equivalent magnitude (e.g. , 
Mather, et al. , 1976b). 

Tlie same cannot be said with confidence, of gravity anomalies, given by 

2AW / AW , A 

Ag = g - T (l + f + in + — - 2fsin^0„ + o {f^Ag}) . (1) 

a \ 2uy } 

where AW is the difference in geopotential between the general point P (ellip- 
soidal coordinates , X on a reference ellipsoid with equatorial radius a and 
flattening f) at the Earth's surface and the regional datum level surface and not 
the gooid. T in Equation (1), is normal gravity computed on the equivalent equi- 
Ijotential ellipsoid rotating with angular velocity , using a value adopted for 
the product of the gravitational constant G and the mass of the Earth M, m being 
given by 


m = (2) 


If a data bank of gravity anomalies is based on IGSN 71, the precision of the 
gravity anomaly field in the context of surface integral evaluations, is influenced 
primarily by the errors in AW, established as a network. It has been shown 
(Mather 1974, p. 102) that the quality of height anomalies estimated from 
gravity anomalies using relations of the form 

= K 


// 


«l)(Ag, AW) da , 


(3) 


1 


where Ap; Is the gravity anomaly and Aw tlio difference of geopotontlal from the 
global datum level surface at the clement of surface area da which Is at a geo- 
centric angular distance from the point of computation, Is a function of the 
wavoleng:th of the errors In the global gravity anomaly data setp Satisfactory re- 
sults are obtained In practice only when the errors In such a data sot decrease 
rapidly as a function of wavolongtlu For example, In tlie case of the gravity 
anomaly data set for Australia - AUSGAD 76, the long wave error sources are 
assessed as being the following (Mather, ot al. , 1976b, p. 79, ot soq,); 

a. Frrors with amplitude 0, 15mGal and wavelength 7000 duo to residual 
errors In the adjusted Australian lovollliig survey. 

b. A constant error of ±0. OOmGal duo to the gravity value adopted for the 
National Base Station at Sydney not being correct, 

c. Errors with amplitude 0,2mGal and wavelength of 70001cm due to resid- 
ual errors In the Australian National Gravity Network, 

A fourth significant source of error when such data banks are used in global 
solutions, Is the effect of the adopted datum J.ovel surface not necessarily coin- 
ciding with the geoid to better than ±lm causing systematic effects of iO, 3 mGal 
in the entire gravity anomaly databank computed on this datum. This statement 
presumes that a world-wide definition has been adopted for the geoid with a res- 
olution to at least ±10 cm. as discussed in (Mather 1977, Sec. 1). 


Tile lack of a global coverage of surface gravity data (e. g. , Bapp 1977, p. 3) and 
tile questionable quality of oceanic coverage which is widely conceded as being up 
to an order of magnitude inferior to land based data, result in the use of combin- 
ation solutions for the geoid. The input data for such solutions are the following; 


a. Satellite determined harmonic coefficients C«nm of fhe gravity field to 
some degree n' (= 20). 

b. Several regional gravity field determinations with representation on, 
say, a lOltm grid - freely available in Europe, North America and 
Australia; providing surface gravity field coverage within 20® of the 
point of computation in selected parts of the regions mentioned. 


The computation can then be carried out using the truncation function technique 
first proposed by Molodensidi (Molodenskii, et al. , 1962, p. 146) using a relation 
of the form 


r - + K' Qn(<^o)Agn + K /*/f(\i/) ‘I> (A g, AW) da , 

n=o /y. 


(4) 


' P ‘^^0 


2 


whore Agn Is Uie n-th dogreo harmonic in Ag, is the outer radius of the cap 
over whicii surface integration Is undertaken, J'r being the residual contribution 
duo to truncating the outer zone effect to representation in terms of a finite sot 
ofhannonicB to degree n*. Several variants of this method are known (e.g. » 
Mather 19C8; Marsh and Chang 1976), 

In the past, the ultimate tost of the accuracy of such solutions haf?.' been obtained 
on Comparison with astro-geodetic determinations. The oomparisuns should bo 
performed between independent estimates of the height above ellipsoid. In prac- 
tice, Stokesian solutions are compared against so-called astro-geodetic geoids 
which can have distortions of up to 2 m In mountainous country. The resulting 
discrepancies over continental extents have root moan square (rms) values of up 
to ±lim, after allowing for datum transformation (e,g, , Mather 1970; Mather 
1975a). Obviously, effects which are linearly variant over the region of compari- 
son are not reflected in the rms residuals. 

Tlie limited precision of astro-geodetic determinations due to the high cost of 
such surveys and the local fluctuations in th; grades of level surfaces, precludes 
any possibility of making more exacting tests on the quality of gravimetric de- 
terminations by those methods. 


2. THE ROLE OF THE GRAVITY FIELD IN OCEAN SURFACE DYNAMICS 

A complete treatment of this problem is given in (Mather 1978a). It can be sum- 
marized as follows? The dynamics of the surface layer of the oceans are defined 
by the differential equations 


- r^:2 = - g 


bL 


9Pa 


Ox, 0X| 


+ Fj + 0 {f Kj} , 


(5) 


and 


x*2 + fxi =-e- — +F2+0 

Oxj Pw 0X2 


|G) 


where f is the Coriolis parameter, given by 

f = 2oi sin0 , 


(7) 


p,i is the atmospheric pressure, (x, , X 2 ), (x, , X 2 ) and (F, , P 2 ) are components 
of accelerations, velocities and frictional forces acting on the surface layer of 
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the oconns along tho Xi and X 2 axes oriented oast and north respectively In the 
local horizon, 

'riio major non-equatorlal quasl-statlonary currents (o. g, , tho Gulf Stream, the 
Kuroshlo) have velocities In excess of 10^ cm S’* maintained by gradients of tho 
quasi -stationary dynamic sea surface topography (SST) defined as the helglit 
of tho sea surface above tho geoid. In mid-latitudes, a steady state current of 
1 cm s’* is maintained by a quasl-statlonary SST gradient of 1. 05 cm per XO^ Ion, 
tho current being deflected in tho direction of the SST contours. The wind ve- 
locities and/or atmospheric pressure gradients needed to maintain currents like 
the Gulf Stream are at least an order of magnitude larger than the strongest 
measured under extreme conditions at the surface of the Earth ns summarized 
in Figures 1 and 2. 

Tho spectrum of SST is illustrated in Figure 3. The quasi-statlonnry component 
in mid-latitude regions where fast flowing steady state currents occur, can be 
expected to bo 4-5 times larger than the time varying constiluonts, tho overall 
magnitude of f 5 not exceeding ±1 '/am, 

Remote sensing techniques provide the only plausible means of synoptically mon- 
itoring the dynamics of the surface layer of the oceans. An Indirect method of 
doing so can be developed from an analysis of infrared Imagery. Tlie primaiy 
difficulty is tho interpretation of relative measurements to provide an absolute 
scale. A second limitation is imposed by cloud cover, Tlie greatest apparent 
strength of Infrared imagery is in tracking eddies. However, tho temperature 
structure Within an eddy is complex and does not lend itself to straightforward 
interpretation (e.g, , Cheney and Riehardson 1976, p. 145), especially in shallow 
seas. 

Tlie radar altimeter provides direct estimates of the position (and hence the 
height) of the instantaneous sea surface and is unaffected by cloud cover, though 
refraction corrections to the measured range may be more uncertain under such 
conditions. If the measurement were reduced to a height (D above a reference 
surface and if the height of the geoid above this same surface were N, it follows 
that 

?s = - N (8) 

The role of the gravity field in ocean dynamics modelling is somewhat different 
from that in tho solution of the geodetic boundary value problem. In the latter 
case, the objective was the geometrical mapping of a surface at which the meas - 
urements were made. In tho ocean dynamics application, the shape of the 
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Figure 1. Equations of Morton at Ocean f^rface 
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Figure 3. Guestiniater of the Space — Time Spectrum of Sea Surface 
Topography Determinable by Satellite Techniques 


l)oundtng surfnco Is lmo^vn, It Is roqulrod to goomotrlcally map (in concept) a 
lovol surfnco In oconn nroas undor olroumstnncoB whore no moasurements have 
boon dtrootly mndo In relation to It. 

Consoquontly, attempts to find a moans for dotonnining SST from a solution of 
the boundary value problem, without making assumptions about the nature of the 
SST, have not been successful (Mather, etal, , 197Ca), In formulating a method 
for dotorminlng ?s from satellite altimetry, it Is not considered desirable to 
assume characteristics for the global distribution v/f ?s processing the data 
prior to obtaining a solution, Consoquontly, the only sources of data on the 
Earth's gravity field which are Independent of any relationship to the geometry 
of the sea surface, are satellite orbital analysis and satelUto-to-satolllte track- 
ing (o, g, , Vonbun, et al, , 1977). No reasons exist at present for assuming that 
gravity field models deduced from a global network of tracking stations to x 
parts in 10** will not achieve a resolution of x parts In 10® when transformed In- 
to long-wuvo components of goold heights through wavelengtlis which are a func- 
tion of: 

• satellite flying height; 

• the nehvork of higher satellites used in satellite-to-satellite tracking; 
and 

• the uvorage distance between stations in the tracking station network, 

as can bo seen from a stuc'/ of crossovers of GEOS-3 data (Mather, et al. , 1978b, 
Table 2). 

Tlie iteration between purely satellite-determined gravity field models like the 
odd-numbered Goddard Earth Models - e, g, , GEM 9 (Lerch, ot al, , 1977) and 
dynamic solutions using a more widespread network of reliable laser ranging 
systems, can be expected to provide improved estimates of the low degree har- 
monics Canin of the global gravity field through some degree n' (equivalent wave- 
length of approximately 10^ km) using global networks of ^tlO cm tracldng sys- 
tems. The possibility exists that the minimum wavelength resolved could be 
decreased to about 500 1cm if data from aatellite-to-satelllte tracking were in- 
cluded in the solution, using data collected during a mission of tlie GRAVSAT 
type, provided the data was of sufficient precision. 

Solution techniques for the recovery of SST from radar altimetry data must re- 
flect the band limited nature of the geoid height signal from satellite determined 
gravity field models. 
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The role of purely yeodotlc technjifiuos In ronnoto sensing ocean dynomtes Is pri- 
marily that of providing values of fjj from the radar altimetry data. The signif- 
icance of this technique lies in Its potential to synoptically monitor changes In 
the shape of the sea sur face, 'Ilie geodetlo results have the potential to define 
a four dimensional fi’ame of referonce for sub-surface ocean dynamics, In 
addition to providing information for modelling the dynamics of the surface 
layer of the oceans. 


3. RISMOTE SENSING fs RADAR ALTIMETRY 

As discussed by Mather, et al, (1976a), it is not possible to use the boundary 
value problem approach in determining SST (i*.;) from sea surface heights f' 
above the reference figure because no data can be unambiguously related to the 
geoid at the desired level of precision without making unwarranted assumptions 
about the magnitude and distribution of the SST, 


Tlie spectrum of ocean surface dynamics lends itself to convenient sub-division 
into two components in the context of remote sensing. The first is quasl- 
stationary in time { ) during the period of data acquisition, Tire time varying 
component fst Is expected to have a magnitude which is about one fifth that of 
?so for reasons given in Figures 1 and 2. Some dominant contributions to ?so 
have hcon recovered from GEOS-3 altimetry and have substantially the same 
magnitude as obtained from oceanographic surveys. The ssmoptic variations in 
the SST are likely to have periods of greater than 2 months through wavelengths 
greater than 10^ km. 


Consequently, the precision required in monitoring the synoptic variations in 
is almost an order of magnitude greater than that needed to establish a gross 
model of qnasi-statii'uary SST. Furthermore, there is no necessity for maps 
of hydrostatically determined quasi-stationary SST (e.g,, Levitus and Dort 
1977, p, 1283) to agree with satellite determined models to better than ±20 cm 
as the former are, in essence, averages overlong periods of time while the 
latter represents neo-synoptic monitoring of the phenomenon. 


The satellite altimetry data from either of the altimeter-equipped spacecraft 
GEOS-3 or SEASAT-A, due for launch in 1978, are in the form of profiles. In 
the former case, the profiles seldom exceed 20 minutes in time due to the ab- 
sence of on-board recording facilities, SEASAT-A will sweep out a 25“ grid 
every day ( Pig. 4). No data is collected outside certain bounding parallels (65“ 
in the case of GEOS-3, 72“ in the case of SEASAT-A), The orbital periods of 
the two satellttes are different, GEOS-3 has an orbital period of lOl, 79 min, 
which results in a daily offset of 1250 km, an n“ x n“ grid being generated every 
25 /n days. The SEASAT-A orbit is planned so that the daily offset is approxi- 
mately 20icm with the groundtracks repeating themselves every 4 months or so. 


OF POOR (JUALllt 




Figure 4. The Ground Track of SEASAT-A With Groumi Station Coverage 





Consequently, as illustrated in Figure 3, SEASAT-A will provide little Informa- 
tion on synov>ttc variations in ocean circulation with wavelengths boLwoon 10^ km 
and 2f)00 Ion and periods loss than 4 months - a potential drawback when attempt- 
ing to model the energy transfer between the variable wind fields and the surface 
layer of the oceans, assuming the former to operate with stmilai Jonstraints in 
space and time. 


4. BASIC RELATIONS 

"fhe following is summarized from (Matherl978a). It is assumed that a harmonic 
representation of the gravity field of adequate precision is available to some de- 
gree n'. It has previously been assumed (e.g. , Mather 1974, p. 90) that the de- 
sired precision in each coefficient C^nm defining the gravity field model was about 
1 part in 10^. Table 1 lists the normalized coefficients fsann, (n < 5) in the sur- 
face splierical harmonic representation of the hydrostatically determined quasl- 
stationary SST based on data confined to the oceans lying between the parallels 
65"S and G!)“N (Mather, et al, , 1978b, Sec, 7). Also listed are errors ecanm in 
the coefficients Caimi of GEM 9 (Lerch, et al, , 1977, p, 52) with their linear 
equivalents. A study of the signal-to-noise ratio (fsanm/ecanm ) shows that con- 
ditions are favorable only for the recovery of the coefficients ‘isiii » ?si 20 » fsi30 
fsHO and, possibly S’sieu • Table 2 lists the root mean sqiiare error per degree in 
GEM 9, which indicates that there is no possibility of recovering any informa- 
tion on SST from present-day gravity field models with wavelengths less than 
10*’ icm. Tlie desired precision in the gravity field model for qiiasi-stationary 
SST determinations is therefore 0, 2 parts in lO’ through wavelengths greater 
than 10^ km, and hopefully, greater than 500 ion if satellite-to-satellite tracking 
methods can provide the necessary precision. 


The desired resolution of the gravity field model for synoptic monitoring SST 
variations, is at least five times more exacting, but only through wavelengths 
which cause perturbations of ±1 cm in the radial component of orbital position, 


The geopolential W exterior to the Earth's atmosphere can be represented in geo- 
centric spherical coordinates (E, 0 , ^ ) by the relation 


W = 


GM 

il 


E i 


n=0 



^anm ^annu ” ^ ^ 


(9) 


where Canm are spherical harmonic coefficients of degree n and order m, S„n„i 
being surface spherical harmonic functions defined by 

^i iitii ~ (sin0) cos mXj ^ 2 nni ~ ^nm (sm0) sin in\ . (19) 
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TaWle 1 


KactorH Influeni’lnn Determinations of (Aiasi-Stationary Dynamic S» a 
Surface Topojjraphv From Satellite Altimetry — 'Phe Sltmnl 
the Noise per Coefficient (Cfanm ) 
and the Sijjnal-to-Noise Ratio (f»anm /©( onm) 











Table 2 

Decree Variances 


Degree 

(n) 

Sea Surface Topography* 
(cm*) 

GKM 9 Error** | 

(KGal cm)* j 

1 

1 

528.6 

” 1 

2 

2170.3 

15 

3 

107.0 

138 ' 

4 

109.7 

83 

6 

23,7 

391 : 

6 

31.1 

266 1 

7 

3.9 

865 

8 

4.5 

619 1 

9 

2.7 

x435 

10 

0.7 

1193 

i ” 

0.4 

2415 

t 

12 

0.1 

1815 

1 

0.1 

2452 

14 

0. 0(4) 

2049 

la 

0. 0(7) 

2575 

16 I 0.0(1) 

2248 1 


♦Based on a constrained solution to (16,16) from 5® x 5® area means In ocean 
areas only fit to data » *9 cm. 


** Based on comparisons with surface jjravitv (Lerch, et al, , 1977), 
Comparisons with altimetry are a factor of two better on the average. 
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The quantity a in this development, is the equatorial radius of the ellipsoid of 
revolution which best fits mean sea level for the epoch of the altimetry. This is 
obtained by analyzing the heights of the sea surface in relation to an adopted 
reference surface in which the ellipsoidal radius (Mather 1974, p, 91, et seq.) 
is ao. The change da in a is obtained from (Appendix, Equation (A-12)) by mini- 
mizing the residuals 


V 


{■' + a df - da - E E fsaim 


al tn > 


m=0 a=I 


( 11 ) 


df being the change ‘in the flattening between the sea surface and the value of f 
obtained from CI 20 « The question of permanent Earth tide effects on C'no is 
discussed in Section 10. 4. The final set of terms in Equation (11) allows for the 
first degree harmonic in the quasi-stationary SST. 

The value of da obtained in practice (e, g, , Mather, et al. , 1978b, Sec, 7) is 
based on a sample which is banded in latitude due to the distribution of the al- 
timetry. The resulting ellipsoid of "best fit" will not be representative of the 
global oceans. This highlights the difficulty of defining the equatorial radius of 
an ellipsoid which "best fits the geoid globally, " There is no way of sampling the 
geoid globally. Land areas could be introduced into the definition by using the 
geopotential differences between tracking stations and the regional mean sea 
level (MSL) datum. The adoption of such a procedure will no longer define an 
ellipsoid which best fits the geoid in ocean regions (Mather 1977, Sec. 1). 

Neither will it provide the additional coverage needed in higher latitudes in the 
short term. 


The approach adopted to date in analyzing GEOS-3 data (Mather, et al. , 1978a, 
See, 6.1) for the selection of a level surface as the geoid is the following. The 
oceans are treated as lying entirely within parallels 65“S and 65°N, All inland 
seas like the Great Lakes and the Caspian are ignored when sampling f' using 
Equation (11). The resulting level surface is defined in terms of the potential 
Wo of the geoid consistent with the value of the product of the gravitational con- 
stant G and the mass of the Earth M included in the adopted gravity field model. 


The geopotential W in Equation (9) cannot be downward continued through the 
atmosphere using a spherical harmonic model. If the atmospheric potential 
is defined at all points exterior to the geocentric sphere of radius Ra enclosing 
the Earth's atmosphere, by a relation of the form 


GM ^ /R^\n " 2 

'^sanm Sanm > 

*'■ n=0 \ *^ / m=0 a=I 


(12) 
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it Is possible to tloClne the potential Wc of the solid Earth and oceans but oxclud- 
inK the atmosphere according to the relation 


Wn 


GM ^ / a\'> " 


= T E i i: i: 


anm ‘^anm 


+ Ur, 


(13) 


11=0 


m=o «= 1 


where U,. if the rotational potential and 






V. 


sail 111 


(14) 


W^; satisfies Laplace’s equation in the space exterior to the Earth's surface 
(air/sea interface at sea). 

[Ills now possible to consider the disturbing potential T" in relation to a higher 
reference model whose potential U in space exterior to the Earth's surface is 
given by (Mather 1974, p, 91, et seq.) 

U = ' j, 2 ^ { 1? ) ^ Uftntji Sournr + Uf . (15) 

nf-o \ / m»o a=i 


Note that tlie potential will contain non-zero first degree harmonics if the 
origin of coordinates (R, 0 , X) remains at the geocenter. It has been shown that 
enforcing the condition 

c'l'io = c;;, = c^'n = 0 (IG) 


will introduce errors less than ±5 cm due to the non-coincidence of coordinate 
origins (Anderson, et al. , 1975, p, 33), The height anomaly f' on this higher 
reference model has a magnitude which is one order smaller ithan die value f in 
relation to a rotating equipotential ellipsoidal reference model, being given In 
ocean areas by the relationship 


GM ^ / a \ « ^ ^ „ 
f f “ T~~ X/ ) X/ X ^■'Ofnin ^anni « 

‘'o' n =2 \‘'o/ m=o a=l 


whe re 


Uoiim " U^nni except when oi = 1, in = 0 and n ~ 2 , 4 and 6. 


(17) 


(18) 
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nnd G|f,o have been correoted for the effect of the ellipsoidal flattening im- 
plied in the value of C '120 by the relations (Appendix, Equations (A-2G)) 


Cl., 


I‘50 




siii'la 


'MO 


1 - 


4 lu sina 


21 q2(a)J 
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and 


C)60 ” C{f,o 


sin"a 


2 m slim’ 


9 ci2(a) 


- I 


( 21 ) 


the relationship between C', 20 , f, ot nnd q 2 (a) being given in the Appendix, 7, 
for all practical purposes, is normal gravity on Uie Implied level ellipsoid Im- 
plied in the higher reference model, is related to the disturbing potential T" 
in relation to tho higher reference model by the relation (Mather 1975, p, 72)s 
Mather 1978a, Equation (19)) 


fpf/ 


- U = 


_ CM 

” I? <0 

n=o 



^ Can,,, 


emm 


1 


( 22 ) 


= (Wo-U«) + - 7{-s - V. 


where W,, is the potential of the geoid which is determinable from satellite al- 
timetry consistent with the higher reference model (Mather, et al. , 1978a, Sec. 
G.l), Ui) is the potential on the surface of the equipotential ellipsoid implied in 
the higher reference system and V is the potential of the atmosphere, 'Ilie first 
expression holds in space exterior to the Earth's surface while the second holds 
at the surface of the Earth. Equations are also given in (Matlier 1978a) for the 
gravity anomaly and solutions of tho geodetic boundary value problem which take 
into account the unknown quasl-stationary SST, the complete zero degree effect 
and conditions of continuity of Hie geopotential in the space exterior to the Earth 
given tho existence of the atmosphere. 

5. PRACTICAL SOLUTIONS 

Numerical results from GEOS-3 data have been reported by Mather, et al. , 
(1978b) where selected low degree harmonics of the quasi-stationaiy SST have 
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boon osHnmlod from a Rlobnl analysis of Iho available GEOS-3 altimetry, llio 
following assumptions are made when using Equation (22) in practice. 

• 'Hie low degree harmonics of the gravity field aro Itnown to some degree 
n" with an error loss than ± ao(= 10’® in the GEOS-3 study), 

'rims 


^ofiim “ ^anm “ ^aniii • (23) 

where and CQ„m dofiend In Equations (14) and (18) respectively. On 
adopting surface harmonic models for i’s (eoefficlonts fsanm )> V (coefficients 
Vj,an„i) and (cooffiolonts ?aii„t), It follows tliat 

fsanm “ fftnm 
El, being the geocentric distance to the sea surface. 

It has been shown (Matlier, et al, , 1978a, Sec, G. 1) that for all practical pur- 
poses, the netoffeotof the terms witliin parentheses should be zero if the value of Ci,,,,, 
used in forming botli and were free from error. In such a case, 

fsanni “ fanm ‘ (25) 


CM 


(jr) I Cm- 


C 1 + 

^annr gamn 


I 


(24) 


Numerical solutions are not so straightforward due to the poor signal -to-noise. 
This problem is dealt with in depth in (Mather, et al. , 1978b, Sec. 7) as are 
the procedures for modelling the SST (ibid. , See. 8). As discussed in Section 
4, it is judged that only five coefficients can be recovered from a perfect de- 
toimiination of T and the GEM 9 gravity field model. Other coefficients can be 
estimated but the level of uncertainty is much greater, being a function of 

• the uncertainty in the value of Uie GEM 9 coefficient; and 

• the magnitude of the coefficient (Tables 1 and 2). 

Table 3 lists provisional estimates for the coefficients fganm from the 

analysis of GEOS-3 altimet er data between parallels 65®S and 65®N and the esti- 
mated errors. For detaiU of the method of analysis, see (ibid.). These values 
are based on the following set of constants: 


1 ? 



Table 3 


Proliniinary Estimrtos of the Dominant Features of Quasi^-Statlonary 
Dynamic Sea Surface Topography from the 1977 GEOS-3 Altimeter 
Data Bank and Goddard Earth Model (GEM 9) 

(For Method ‘Of Selection, See Table 1) 


Coefficient 

Normalized Value (cm)* 

Oceanographic 
From Sterlc 
Anomalies 

Geodetic 
I'rom Altimetry 
and GEM 9 

110 

4«.9 

-150 ± 15** 

111 

-21.8 

13 -M 15** 

211 

2.4 

20 -b 16** 

120 

-46.2 

-43 ±6 

130 

6.7 

+7 ± 10 

140 

-9.5 

-18 ± 15 

160 

4. 4 

1 ± 15 


+Based on an analysis banded between 65“S and 65“N 
**GEM 9 does not include a first degree harmonic, system origin being coin- 
cident with geocenter. These numbers also include the displacement of the 
estimated location of the geocenter in GEM 9 tracking station coordinates 
from the center of the ellipsoid which best fits the sea surface. 
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c « 2,997 924 58 X IO>“cms-‘ 
CM - 3,986 004 7 x 10^0 cin^ sr^ ^ 


( 20 ) 


micl Oio f?oomotrical cnlibrntion of the GEOS-3 altimeter performed by Martin 
and Butler (1977), using' the Equinox Data Sot (Mathex*, et. al, , 1978b, Sec, 4), 

'llie geopotontial of tlio geoid for the epoch 1970. 0 from the GEOS-3 altimetry is 

1 

^ Wo = 6,263,283,8 ± 0.4 kGal m. (27) 


'rho ellipsoid of revolution which best fits llio sea surface is defined by the 
parameters 

II = 6,378, 1 39,9 ± 0.4 m: f= 1/298.237 ± 0,003, (28) 

the latter being consistent with the value of -43 cm for the normalized value 
f.si 2 n* These provisional values are subject to revision with Improvements 
in the orbital ephemeris of GEOS-3, The contribution of the permanent 
Earth tide has been eliminated as described in Section 10,4, 


fi. FUTURE TO ENDS 

(a) G cavity Model Improvement 

I'lirther progress in determining quasl-stationary SST is dependent on the reso- 
lution with which the global gravity field model is defined. This information 
must, In the final analysis, come from orbital analysis and satellite-to-satollite 
tracking as envisaged in Section 3. The greatest limiting factor at the present 
time is the absence of a global network of high precision tracking systems. In 
the short term, the satellite altimetry can play a role in gravity model 
improvoment, 

As some dominant harmonics of the quasi-statlonary SST have been evaluated, 
it is possible to remove Uielr effect from the data base of t' obtained from 
GEOS-3 altimetry between 65“S and 65‘’N, assuming that sufficient passes of 
altimetry are available to eliminate the effect of the tides, as well as mesoscale 
variations in the SST on the data, Tlie objective is die evaluation of the coef- 
ficients of a gravity field model with a resolution of 1 part in 10® (i, e. * GkGal 
cm) in each harmonic. ITie best available model at the present lime (GEM 9) 
appears to satisfy this precision requirement to degree 4. Thus, holding these 


ImnnonlcB heavily constrained. It is possible to attempt a gravity model Im- 
provement using observation equations of the form {Mathori ot al. , 197Gb, p, 

42) 

GM i / a A ^ 

2 J (n ) zJ iC ^aiim Sa«m “ 'I'J'sr “ ^ > (29) 

*'0 ii“2 \‘'p/ m“0 oti*\ 

on the GEOS-3 altimotry data bank, whore 

li II 2 

fsr “ ^ ^ 2 ?sanm ®anm "•* “ ^o)“ YJ Z')' > (20) 

n“0 m»o a“i 

whore n ropresonts the maximum degree to which the coefficients of the SST arc 
kno^vn to -16 cm. 

Ii\irthor progress is possible only if it is assumed that the magnitude of higher 
degree coefficients of the quasi-stationary SST are significantly smaller than 
the errors in GEM 9. Tliis refers to the set of coefficients Canm comprising 
the current gravity field model which has to be updated and is used in defining 
the higher reference model (Equation (15)), Table 2 sets out the error degree 
variances in GEM 9 and the degree variances In the oceanographically deter- 
mined SST to degree 16, showing that considerable scope exists for refining 
gravity field models from altimetry data corrected for the effect of the dominant 
features of the SST. 

Assuming that the altimetiy data banlt were of adequate precision, it follows that 
the resulting gravity field model should also have a resolution of 1 part in 10® 
if there were no aliasing effects due to the data being restricted to ocean areas 
only (33,902 out of 64, 800 oqul-angular 1* x 1" squares) and banded in latitude 
between OG^S and CG^N, This aliasing problem is discussed at length in (Mather, 
et al. , 1978b, Sec. 8), 

Tlie goal of complete areal representation can be improved by forming observa- 
tion equations using gravity anomalies for land areas in the relation (ibid. , p. 

42) 

GM /iX** A A 

2 ^ Canm - Agci = V, (31) 

^ 11=2 X*'/ ni“0 a=i 

where the spherical harmonic series is evaluated at the surface of the Earth, 
the value of observed gravity (g) is used to compute Agj, defined by 


20 


( 32 ) 


Aga ° Ag - fiT + 6g., - - I (Wo - Uo) 

n n 2 

’iT Xi/ S fsanin ^anm ■** o{n*/K) , 

' iinfl iii“0 aol 

Ag bBlnt? dofinod by Equation (1), By being the oorreotion to Ag v'hen using the 
higher roforonce model (Mather 1974, p, 96), Sg^ tlio correction for the atmos- 
phere (Anderson, ot al. , 1976, p. 26) and is tlie deflection ol the vertical. 

It is assumed that the dominant harmonics of the quasl-statlonary SST as do- 
tormlnod earlier, adequately model the h-^ight of MSL at the regional lovolllng 
datum, For some studies of this problem, see (Mather, ot al, , 1978a, Sec. 7), 

It is doubtful whether gravity field models developed in this manner, can play 
any role In determining higher degree harmonics In the global quasl-statlonary 
SST foi' obvious reasons. However, the use of such n field as the input model In 
the. ro-nnalysls of orbital data for a ftirther reflnomont of the gravity field model 
Is the next oljvlous stage In the proeess of improving the gravity field to the pre- 
cision of 2 parts in 10*^ required for ocean dynamic modelling, 

(b) Determination of Shortwave Contributions to the Qunsi-Stationary SST 

'I’hose short wave contributions are defined as those with wavelongUi less than 1!, 
whore C is the wavelength of the highest full harmonic of die gravity field which 
porturljs the orbits of near Earth satellites in excess of the noise level of the 
tracking. For a discussion of this subject, see (Mather 1978a, Sec. 7). This 
discussion proposed a technique based on the solution of the Inverse of the geo- 
detic boundary value problem, formulated to take into account the departures 
of the sea surface from the geoid (considered known with wavelengths greater 
than il), 

llie equations comprising this solution were the following, 

TN'.' F « 2 12 

Age + - (W,j - Uo)-T fsanm ^anm 2-/ ^saim ^ai 

L 11=2 111=0 a®! m°o a=i 

(33) 

- 0 
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whoi'o 


f\f') ° n(2n + 1) I’„o (cost//) , (34) 

I1“2 

^ being the goocontdo angular distance of Uio element of surface area d a from 
the point of computation P, 

Nc “ f' “-^ (V-ST"). (35) 

Age = Ag - 5? + 5ga - ~ + 5Ag" , (3G) 

6t", 6Ag" being the changes in T" (Equation (22)) and Ag" defined by 
cr ■ 53 ('jr) ^anni ^«nm • 

no 2 \‘'/ ni“0 tto 1 

in Equation (31), between the Earth’s surface and the minimum geocentric sphere 
onclofiing llie Earth’s topography (the Brlllouln sphere) of radius R. 

'rhe principal difficulty In using this equation In practice Is the requirement for 
appropriate resolution in Ag (2 parts in 10^) through wavelengths of Interest, If 
It Is to play a role in ocean dynamic modelling. The data Is required within 
approximately 500 km of the point of computation, with values of required 
everywhere, including land areas. It appears optimistic to expect results from 
this technique in the foreseeable future as surface gravity data of adequate pre- 
cision will not be available for the task. Thus surface gravity Information does 
not provide a viable basis for recovering s’'ort wave information on fs. 

An altoinato means, in principle, for obtaining improved definition of the gravity 
field through wavelengths between 5001cm and 20001cm with a precision adequate 
for ocean dynamic modelling, is by satellite- to-satelli to tracking of a low-flying 
satellite with appropriate precision from a network of spacecraft In synchronous 
orbit. Research in this area is still at a vei*y early stage. ITio desired resolu- 
tion is better than 10"^ cm s*2 with wavelengths in the range mentioned above. 

Tills Is almost two orders of magnitude smaller than the resolution achieved to 
date ( Vonbun, et al. , 1977*, Marsh, etal, , 1978). More recent analysis appears 
to indicate that it is possible to recover accelerations at the low flying satellite 
altitudn to J;2 x lC'^cm s’^ from range-rate data with a noise level of J:0, 2mms'' 
(Marsh 1978). Tliese results are based on minimal quantities of data and it 
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would bfl pessimistic to a<jsurao that they cannot bo improved upon substantially 
by lucraglng over long periods of time for the c’ sired effects of intermediate 
wavelength* 


7. SYNOPTIC MODELLING 

Synoptic modelling of tlio circulation patterns of the surface layer of the global 
oceans on a dally basis with wavelengths greater than 2600km is a very real 
possibility with the launch of SEASAT-A. The computer requirements are well 
within the capability of the Goddard network. The principal limitation at the 
present time is the precision of radial orbital determination using the present 
tracking network as well as the weakness of current gravity field models through 
harmonics other than low degree '/onals. It is not beyond the realm of possiijility 
that the altimetry data available from the GEOS-3 mission, appropriately cor- 
rected for the dominant features of the quasl-stationary SST as determined in 
Section 6, can bo combined with current models and tracldng data to produce a 
gravity field with a resolution of 1 part in 10*^, 

'Iho pi’incipal difficulty is the aliasing effect produced by the altimetry data being 
confined to the region between parallels 65®S and 65“N and the lack of data of 
ecfuivalcnt precision on land. Another factor which will help resolve this problem 
is the increase in the density of tracldng stations in the high precision global 
tracking network with the passage of time. 

The altimetry data can also be used with local high precision tracking data to 
study regional circulation patterns with periods greater than 1 month. Studies 
of eddies and other features with vertical magnitude in excess of 20 cm. have been 
reported from GEOS-3 altimetry (Mather and Coleman 1977; Leitao, et al. , 1977; 
Mather, et al. , 1978c). All these studies were done without benefit of precise 
orbit determination and lose some part of the spectrum of the SST in space and 
Lime as a consequence. Tliese evaluations are based on determinations of changes 
in Urn SST between epochs (t = t) and (r = t + dt), using relations of the 
form 

65-s = rCt + dt) - r'(t), (38) 

where is established directly from the altimetry and the orbits together with 
the enforcement of crossover constraints to eliminate any unmodelled force 
field errors in the orbit integration. 

It can be stated with confidence that the use of geodetic techniques on radar al- 
timeter data Is likely to produce information which Avill make possible tlie syn- 
optic monitoring of the dynamics of the surface layer of the oceans. 
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10. APPENDIX 


10.1 Relations Between Geodetic and Geocentric Pax’ameters to 1 Part in 10 


See Figure A-1. The geocentric latitude, <f > , the geodetic latitute 0g , the geo- 
centric distance R to tlxe ellipsoid and the radius of curvature iv in the prime 
vertical, given by 

u = a/(l - c2sln^0g)''^ , (A-1) 


26 




Hg'uro A-1. The Meridian Ellipse 


where Ihe meridiun ellipse is defined by the equatorial radius a and eceentrieity 
e, by the relations 

Kmii0 » ('sin^g - I'c^sin^g = i^(\ -O^sin^g, (A-2) 


and 


R cos (t> = V cos 0g . 

the Qaltening f being related to e by the relation 

c*2 = 21(1 - f/2) . 


Thus 


tan 0 = ( I - f)^ tan 0g . 


(A-2) 

(A-4) 

(A-5) 
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If 


6(^ = 0g - (^ = o {f} , 


(A-6) 


it follows that 


Ian 8<f> ~ 8(j) -i- o 



tan <^g - tan^ 

1 - tail 0g tan 0 


(A-7) 


On using Equation (A-5) in Equation (A-7) and expansion using the binomial 
theorem, it follows that 

60 = r sin 20(1 + f(3/2-2sin2 01 + o (f3) . (A-8) 


On using Taylor’s theorem, the following results are obtained; 


sin20g = (I +4f+ lOt2)sln2 0 - 4f(l + 13/2f)sin^0 + I6f2sin’’0 + o {f^} , (A-9) 


and 


sin‘^0g= (1 +8f)sin‘^0 - 8fsin^0 + o . (A-10) 


The geocentric distance to a point on the ellipsoid at latitude 0g Is obtained by 
using Equations (A-1) and (A-4) in the combination of Equations (A-2) and (A-3) 
to eliminate 0, when manipulation using the binomial theorem gT,es 


R = a 


fO 



+ 2f-)sin2 0g 


r~ 

(5 - nosin'* 0„ 

2 ^ 


13 

T 


f3sin^0g + 0 {r'*} 
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(A-11) 


wliich is the expression for R in geodetic coordinates. 

The expression in geocentric coordinates is obtained by substituting Equations 
(A-9) and (A-10) in (A-11), when the relation 
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(A-12) 


K = a 1^1 - f ^1 + I f + 2f2^ sin2 0 
+ -^ f2 -i- ~ sill‘d <p sin^ 0 + 0 1 1'* I 

is obtained, 

10. 2 Tlio Sphei'ioal Harmotilc Kepresontation of the Exterior Potontlal Due to a 
Rotating Eqiiipotontiai Ellipsoid to Ordei* £4 

If a wore tho parameter defined by 

coscx = 1 - f (A-13) 


mid 

sin a - c = (2f - (A-14) 


and if tho equatorial radius a and polar radius b of the equipotentiai ellipsoid 
wore defined by 


a = ccosQcao! b = ccotao,, (A-15) 

where etj, is the value of ot on the equipotentiai ellipsoid, the exterior potential 
of the rotating equipotentiai ellipsoid at the general point (a , u) is given by 
{e.g. , Mather 1971, p. 83) 

CM a2oj2 

U = g 4- q^(a) P^oCsin u) -t Up, (A-lC) 

c 3q2(ao) 


M boing the mass of the ellipsoid, u the reduced latitude (Hgure A-1), Ur the 
rotationiil potentiiil and q 2 (O') is defined by (ibid., pp. 78-9) 

q2(oc) = “ [o.'(3cot^a' + I)-3 cq 1 q:J 
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Dug to rotational symniotry, the potential in sphoi’ical harmonics duo to tliis samo 
rotational oqulpotentliil ellipsoid is also given by 




P2no(sin0)C2„o 


+ Ur 


(A-18) 


Expressions (A-16) and (A-18) are identical, The coefficients Cjno be eval- 
uated by compailng the two expressions at the pole of the i-otating equipotentlal 
ellipsoid when 

R = b, u = f = 90“ (A-19) 

and ali zonal harmonics talce the value 1. In such a cirqumstoiice, 

c c 

tan (X = ~ (A-20) 

b R ' 

from Equation (15). From Equation (2), 

a2to2 = CM in /a , (A-21) 

Further 

« = Liur*(tana) = tana + 2^ (-1)” — (A-22) 

n«i 2n+I 


Using Equations (A- 17), (A-21) and (A-22) in Equation (A-16), simplification 
after re-arrangement of terms, gives 


GM 

^ tan^'ia 

" 

in sin a 2n 

““ 

U = 

1 + zJ (-1)" ■ 

1 

_ ^ — , 

p 2 o(sin u) 

b 

n=i 2n + l 


3qj (a) 2n + 3_ 



+ Ur.(A-23) 


Equating coefficients of the same degree in E after imposing the conditions at 
(A-19) and (A-20) on Equations (A-18) and (A-23) gives 



ain 


c2n 

= (- 1 )" 

2n+ 1 


ni sin 2n 
SqjCa) 2n + 3 


(A-24) 
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whore tho subscript o has been suppressed in the value of a wliich refers to the 
equipotentlal ellipsoid. As Equation (A-15) defines 

c 

•“=smo: (A“25) 


it follows that 



2n + 1 


1 


msina 2n 
Sqjfa) 2n + 3 


(A-26) 


10. 3 The Flattening f and Cjo 


The flattening f is obtained from the second degree zonal harmonic in the spher' 
leal harmonic sexies by the relations 


or 



2m sina 


LiSqjCoi) 
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(A-27) 




”2m sin^ a 
lSq2(a) 


*■ 

■* 3 C20 


'/i 


(A-28) 


where q, («) is defined by Equation (A-17). TMs obvloux' has to bo solved by 
iteration with tlic procedure rather unstable unless sufficient significant digits 
are carried in computations. 


In most computations, a trial value (C2ot) is available for C20, corresponding to 
a value f( for f, A solution procedure is required for computing df from dC2o , 
given ixy 

CIC20 C20 - €201 , (A-29) 


using an expression for clf evaluated in terms of a rapidly convergent series. 

The use of Equation (A-17) enables the complex term in Equation (A-28) to be 
written as 

2msin^a 

m cosset 

15qo(a) 


+ — tan-a + tairce + o 
7 49 


(A-30) 
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Chfuigos d« in a are produced from changes dC 2 o In C 20 according to 
da 


/ 6 1 ' 

2 sinacosa + 3 m cos^asina 1 + — tan^a + taiV'a 

V 7 49 


12 . 

- — sina 
7 


(A-31) 


+ 3 dCao = 0 


wMch is obtiiined by differentiating Equation (A-27) after re-arrangement of 
terms. Then, as 



(A-32) 

(A-33) 


and 


f = ft + df. 


(A-34) 


10. 4 The Effect of tlm Permanent Earth Tide on Ocean Dynamic Modelling 

The coefficients of satellite determined gravity field models Cn^ referred to in 
Section 4 are computed after modelling the effect of the Earth tide. The latter 
has a permanent constitimnt which causes the second degree zonal harmonic of 
the true gravity field wliich influences the quasi-stationary component of ocean 
circulation to be different from Cj 2 o« Yhis problem is discussed in (Mather 
1978b). The incorrect pre-processing of the altimeter data for the effect of the 
Earth tide can cause errors of 13 percent in the value of ?si 20 as determined 
from Equation (24^. 

The magnitude of the error dfj .120 is described by the equation 

J^si20 (^) + 1^2 - h2)P2o(siii5i I (A-35) 


where (R| , R 2 ) s-re tlie distances between the centre of mass of the Earth and 
those of the sun and moon (masses M) , M 21 declinations fij , 62 )> ^21 1^2 being 
the second degree Love numbers, M ( } referring to tlie mean value over the 
period of altimeter data acquisition. 
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Tills cori'oction is only neoossary when sou surface heights from satellite altim~ 
etry httve not been corrected dymimically for the effect of the Earth tide, The 
dynamic correction allows for the change in the instantaneous Earth space po- 
sition of the geoid, defined as the level surface (W = Wq ) due to the influence of 
the Eartli tide. It is essential that the correction made by dynamic rather than 
geometric (i.e. , baaed only on the Love number ha ) when using altimetry data 
in ocean dynamic modelling. 

This will onsuro that the satellite altimetry data will bo referred to the instan- 
ttmoous datum level sui’faco rather tluui some artificial static geoid and thereby 
ensure that the aliasing effect of the permanent tide is satisfactorily removed in 
determinations of SST. The value obtained for is -IG.I cm. 
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